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Abstract. We describe a simple but effective iterative procedure specifically designed to destripe Q and U Stokes parameter 
data as those collected by the SPOrt experiment onboard the International Space Station (ISS). The method is general enough 
to be useful for other experiments, both in polarization and total intensity. The only requirement for the algorithm to work 
properly is that the receiver knee frequency must be lower than the signal modulation frequency, corresponding in our case to 
the ISS orbit period. Detailed performances of the technique are presented in the context of the SPOrt experiment, both in terms 
of added rms noise and residual correlated noise. 
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1. introduction 

Low frequency noise is known to affect all radiometers and 
to induce correlations among successive samples of the mea- 
sured signal, leading to typical striping effects when produc- 
ing sky maps. It is characterised by a power-law spectrum 
S{ f) oc (l/ff, with /3 roughly in the range 1 2.5 depending 
on the noise source. The total instrumental noise power spec- 
trum is usually specified in terms of the knee frequency, /k, at 
which the white and low frequency components of the noise 
are equal: 



-'7 



(1) 



where cr^ represents the white noise variance per second. As 
shown by Janssen et al. (1996), when data are taken from spin- 
ning spacecrafts, where the spin acts as an ideal switch or chop- 
per, most of the low frequency noise can be removed by soft- 
ware, provided is lower than the satellite spin frequency, f^. 
For present and future space experiments primarily dedicated 
to Cosmic Microwave Background (CMB) anysotropy mea- 
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surements, like MAP' and PLANCK^, destriping techniques 
have already been set up to clean the Time Ordered Data (TOD) 
from low frequency noise. The impact of different scanning 
strategies on the quality of final maps has been studied as well 
(Delabrouille 1998; Maino et al. 1999). 

The problem of removing 1 // noise from the data stream is 
even more important for experiments attempting measurements 
of CMB polarization (CMBP), where even a tiny amount of 
residual correlated noise, at a level negligible for investiga- 
tions of CMB temperature anysotropy (AT/T ^ 10"^), might 
distort the signal characteristics: in the most favourable case 
the expected polarization level is only in the order of 10% of 
the temperature fluctuations. The first algorithm specifically 
designed for destriping polarization data has been proposed 
by Revenu et al. (2000), who extended to polarization data 
ideas previously studied to destripe PLANCK anisotropy data 
(Delabrouille 1998). 

The SPOrt^* experiment onboard the ISS (Cortiglioni et al. 
2002; Carretti et al. 2003) is equipped with 4 correlation po- 
larimeters directly measuring the Q and U Stokes parame- 



' MAP homepage: http://map.gs fc.nasa.gov| 

^ PLANCK homepage: http ://as tro. estec . esa. nl/Planck| 

^ SPOrt homepage: http://sport.bo.iasf.cnr.it. 
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Fig. 1. Top: The sky, in Celestial coordinates, as scanned by 
SPOrt in few orbits. Bottom: Pixel observing time (seconds) for 
388 days of data taking. The pixel size is about 7° (HEALPix 
parameter A^i/t/e = 8). 



ters in the 22-90 GHz frequency range, with FWHM = 7°. 
Besides providing polarization maps of Galactic synchrotron, 
SPOrt will try a first detection of CMBP on large angular scales 
({ < 25), where the CMB E-mode power spectrum is particu- 
larly sensitive to the value of the optical depth of the reionized 
medium (Seljak 1997; Zaldarriaga et al. 1997). Such informa- 
tion cannot be extracted from CMB temperature anysotropy 
data, nor from CMBP data taken by balloon and ground-based 
experiments looking at small sky patches. This makes CMBP 
measurements on large angular scale of utmost importance. 

The algorithm presented here is particularly suited to 
SPOrt, acting directely on Q and U Stokes parameter data. 

The SPOrt radiometers are based on InP HEMT ampli- 
fiers, for which the low frequency noise is known to have a 
l//-like spectrum (J3 ^ I) dominated by transistor gain fluc- 
tuations (Wollack 1995), with knee frequencies in the range 
100-1000 Hz. However, when using correlation techniques, the 
knee frequency is reduced by a factor (roffset/rsys)^, with Toffset 
and Tsys the instrumental off'set and the system noise tempera- 
ture, respectively (Wollack & Pospiesalski 1998; Carretti et al. 
2001). The present state of the instrument already guarantees 
a value /k < 2.5 x 10"^ Hz (Carretti et al. 2003), close to the 
goal knee frequency of one tenth of the ISS orbit frequency, i.e. 



The SPOrt scanning strategy is bound to the ISS motion. The 
ISS orbit is tilted by 51.6° with respect to the Celestial equa- 
tor and is characterized by a frequency /□ = 1.8 x 10"^ Hz. 
The SPOrt antennae, looking at the zenit, will cover about 80% 
of the sky under precessing circles intersecting each other. The 
nominal sampling rate is 1 Hz, and a sky map will be provided 
every 70 days in a mission lasting at least 18 months. Further 
details about the SPOrt sky coverage are shown in Fig.^ 
The iterative algorithm described in this paper has been specif- 
ically designed to destripe Q and U Stokes parameter data. It 
is general enough to be valid for any experiment where data 
are taken by scanning the sky at a frequency higher than the re- 
ceiver knee frequency, and where different scans intersect each 
:Other in a sufficient number of points. Measurements of scalar 
quantities are easily handled. 

The description of the algorithm is provided in Sect.|2 whereas 
Sect. 13 presents various tests of its performances and Sect. 0] 
summarises our conclusions. 



2. Destriping method 

The only assumptions our algorithm relies upon is that the ra- 
diometer is stable during the signal modulation period (the time 
needed to complete one orbit in case of SPOrt), so that the in- 
strumental offset does not change significantly during this time, 
and there is enough overlap between different orbits. In such a 
case the noise can be split in two parts: for time scales shorter 
than the orbit period (high frequencies) it is essentially white, 
whereas for longer timescales (low frequencies) it also contains 
the 1 // component, which we model as a different constant off- 
set for each orbit. A simple iterative procedure is then applied 
to remove these offsets from the TOD before map-making. The 
map-making itself consists in averaging the cleaned TOD cor- 
responding to the same sky pixel after rotating them into a fixed 
reference frame. Parallel transport to the center of each pixel is 
also implemented, as described by Bruscoli et al. (2002), be- 
fore averaging, and can optionally be switched on, though it is 
not strictly necessary for the portion of sky covered by SPOrt. 
It would be mandatory if polar caps were included in the ac- 
cessible region. The time needed to simulate 1 year of SPOrt 
data taking with the nominal sampling rate of 1 Hz, destripe the 
TOD and produce Q and U clean maps is about 10 minutes on 
a 1.3 GHz Pentium III equipped with a 4 GB RAM when us- 
ing 10 destriping cycles, and increases by 2.5 minutes for each 
additional bunch of 10 iteration cycles. 

2.1. Formalism 

Each datum consists of a {Q,U) pair as measured in the po- 
larimeter reference frame: the same sky point, when seen from 
different orbits, shows different iQ,U) values even in absence 
of noise. The standard longitude-latitude fixed reference frame 
(Berkhuijsen 1975) is used when considering the real sky emis- 
sion, which is projected into an A^-pixel map. Our M TOD, 
collected during R orbits about the Earth, are used to build an 
M- vector Y of measured {Q,U) pairs: 



£f™' = 1.8 X 10-5 Hz. 



Y X + B OFF + N. 



(2) 
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X is an A^- vector whose elements are the (Q,U) pairs of the true 
sky emission; A is an M x pointing matrix connecting sky 
pixels to observations (Wright 1996; Tegmark 1997a): its ele- 
ments are 2x2 matrices rotating the Stokes parameters from the 
fixed to the polarimeter reference frame; OFF is an /^-vector 
containing a pair of radiometric offsets for each orbit; B is an 
Mx/? pointing matrix connecting orbits to observations, and 
is an M-vector containing pairs of white noise. In more details, 
the elements of the above mentioned vectors and matrices are: 

Q'] 



(3) 
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IS observed 
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(4) 



where a, is the angle between the polarimeter and the fixed 
reference frames at time t. 
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(5) 



(6) 



(7) 



(8) 



2.2. Iteration cycle 

If we knew the sky emission vector X we could use it to obtain 
the offset vector OFF by inverting Eq. (|2}. Of course, estimat- 
ing the noise offsets would not be worthwhile if we already 
knew the true sky emission. However, to start we can guess 
a 0-order sky emission X'"*, the most neutral choice being a 
null vector; then obtain a 0-order offset estimate, and subtract 
it from our measurements to get a 0-order cleaned TOD; from 
this compute X'^*^ the P'-order estimate of the real sky emis- 
sion, and iterate the cycle. 

Starting from Eq. (|2jl and following Tegmark (1997b), the best 
estimators for the offset vector and the signal matrix can be 
written: 



OFF 

X 



Wp„(y - 



■ A - Z) 
B OFF) 



(9) 
(10) 



where, in case of uniform white-noise variance, as expected for 
the SPOrt receivers. 



1 

obs 

A^-' ■ 

pix 



B' 



(11) 
(12) 



with A^obs = B^B and A^pix = A^A diagonal matrices contain- 
ing the number of observations in each orbit and in each pixel, 
respectively. In such a case Eq (10) simply corresponds to ro- 
tating into the fixed reference frame the offset-subtracted ob- 
servations and averaging those corresponding to the same sky 
pixel. When starting from a null guess map, i.e. Z'"* = 0, the 
first offset estimate is a simple average of the measurements 
taken over each orbit. 

In the more general case of non-uniform white noise, i.e. 
when the noise covariance matrix C - {NN"^) has unequal di- 
agonal elements, the matrices Wobs and Wpix become (Tegmark 
1997b) 



Wobs = {B^C ^B)^B^C ^ 

Wpix = {a'c-'a)-'a'c-' 



(13) 
(14) 



meaning that the simple averages mentioned above should be 
replaced with noise-weighted averages. 

In any case the first offset estimate writes: 



The 0-order cleaned TOD are: 
y(0> ^Y -B OFF^°^ 

and the P'-order estimate of the true sky emission is 

Z(') = Wp,x-F<'" 

The cycle writes: 

OFF'-''' = W„UY - A ■ X''^) 
Z^'+i) = Wpix(F - B ■ OFF''^) 



(15) 

(16) 

(17) 

(18) 
(19) 



and must be iterated until the difference between the sky emis- 
sion estimates of order / and / + 1 is negligible when compared 
to the white noise level. 

The proof that the procedure is effective has been ob- 
tained numerically, as described in detail in the next section. 
Destriping measurements of total intensity (or, in general, any 
other scalar quantity) can be simply performed by replacing 
pairs with scalars in the previously defined vectors, and rota- 
tion matrices with scalar 1 . The average sky signal is lost when 
destriping maps of scalar quantities. However, for maps of po- 
larization data like Q and U it can be kept provided the po- 
larimeter reference frame rotates about the fixed frame while 
running along each orbit, as in the case of SPOrt. This is a nice 
feature, especially when measuring foreground contributions. 

3. Simulations and tests 

We test the method with numerical simulations of the data 
stream expected from about one year of SPOrt data taking. 
As real sky emission we consider both CMB and synchrotron 
radiation, both convolved with a 7° Gaussian beam. Other fore- 
grounds are expected to be less important at the frequencies and 
angular scales covered by SPOrt (Tegmark et al. 2000; Lazarian 
& Prunet 2002) and are not included. 

CMB emission is generated by the CMBFAST"* package ac- 
cording to a ACDM cosmological model with Q„, - 0.3, 

* http://www.physics.nyu.edu/matiasz/CMBFAST/cmbfast.html 
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= 0.7, T = 0.2, //(, = 65 Kms-'Mpc-i, Q.b = 0.05 and 
no contributions from gravitational waves. For the synchrotron 
radiation we use Q and U map templates at 22 GHz developed 
by Bernardi et al. (2003) and featuring a P,^, = y/W^W) - 
35 juK. The polarized intensity peak emission is about 130 ;uK. 
The main ingredients of the template are total intensity low fre- 
quency data (Haslam et al. 1982; Reich 1982), from which 
the syncrothron polarized intensity is extracted, and optical 
starlight data (Heiles 2000), from which polarization angles are 
derived. The template frequency has been chosen to correspond 
to the SPOrt channel where synchrotron emission is expected 
to be most important. 

Noise streams to be added to the sky signal are built in the 
frequency space, where their Fourier coefficients are generated 
according to the spectrum defined in Eq. O, with /3 - I, and 
transported to the time domain via FFT. Our code allows us to 
calculate the angle between our polarimeters and the longitude- 
latitude fixed reference frame for each data sample, and to set 
the sampling frequency from the nominal rate (1 Hz) to any 
other rate while keeping the SPOrt nominal istantaneous sensi- 
tivity of 1 mKs'^^ (Carretti et al. 2003). In the tests described 
here the sampling rate was set to 1 Hz, and the SPOrt observ- 
ing time was about 1 year. The final output of our simulations 
are TOD and noise arrays, clean Q and U signal maps to in- 
vestigate the effects of SPOrt partial sky survey, Q and U maps 
containing both signal and noise before and after destriping as 
well as a hit map to evaluate the time spent over each map pixel. 
We use the HEALPix^ pixelisation scheme (Gorski et al. 1998) 
with Nside - 64 to ensure the effects of pixelisation are negli- 
gible when compared to the antenna beam (the pixel size of the 
output maps is ^ 1°). 

We test the destriping efficiency on simulated maps by evalu- 
ating the residual low frequency noise in terms of added rms 
noise with respect to the white noise case, angular correlation 
functions C^(0) {X = Q, U) and power spectra C] (F = E, B). 

3.1. Correlation functions and power spectra 

To avoid edge effects due to partial sky coverage and to effi- 
ciently take into account the presence of noise, even in case of 
non uniform sky coverage, we primarly test our destriping tech- 
nique by studying the two-point correlation functions C^{6) 
{X = Q, U) measured from our simulated maps. We sample 
the correlation functions with about one degree of angular res- 
olution according to the expression: 



c^{e)- ' 



(20) 



Here Af is the content of pixel / of map X, and Wij is the weight 
of the ij pixel pair 



IWiWj 
Wi + Wj 



(21) 



where each pixel weight w, is proportional to the pixel observ- 
ing time and Hij"!' w, = Npix- In the computation the Stokes 

^ http://www.eso.org/science/healpix 
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Fig. 2. power spectrum of our selected CMB model, con- 
volved with the SPOrt nominal FWHM of 1° (solid), compared 
to the average of 100 power spectra measured by either our 
method (top) or ANAFAST (bottom) on simulated maps featur- 
ing the SPOrt cuts. Measured 1-cr bands due to cosmic variance 
are shown as well. 

parameters are properly rotated so that Q is defined along the 
great circle connecting the points in each pair (Zaldarriaga 
1998). The exact value of the angle 6 at which the correlation 
functions are sampled is calculated as the weighted average of 
the angles between pixels contributing to the same measured 
point. As an example, the measured two-point correlation func- 
tions of Q,U maps containing only white noise are on average 
zero everywhere but at 0-distance, where they are checked to 
be equivalent to the average pixel variance. 
The C^'^ power spectra can be recovered from the mea- 
sured correlation functions by integration (Kamionkowski et 
al. 1997; Zaldarriaga 1998): 

cf= r [c^{0)F\^{ff) + c"{0)F]^{ey\i,m.ede 

Jo 
Jo 

where the functions F'^^i^^) have been defined by Zaldarriaga 
(1998). We perform the integration by Gauss quadrature, 
where the correlation function values at the needed points 6 are 
obtained by interpolation, with a second order polynomial, of 
the three closest measured points. This method of measuring 
power spectra is similar to that used by Szapudi et al. (2001), 
though here it has been extended to polarization data. It has 
been tested on full-sky maps of simulated CMB emission, con- 
volved with a Gaussian filter with FWHM=7°, by comparing 
the measured power spectra to those obtained by using the 
ANAFAST package: we get the same results, within numerical 
precision, up to ^ = 25, i.e. in the full multipole range covered 
by SPOrt. On partial sky maps this method is found to work 
better than ANAFAST, thanks to the lack of edge effects: our 
average of power spectra extracted from CMB maps without 
noise and featuring the SPOrt sky coverage matches the input 



(22) 



(23) 
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theoretical spectrum very well, whereas that measured with 
ANAFAST is flatter, as shown in Fig. |2] Furthermore, when 
trying to extract the real sky signal from noisy maps this 
method allows for pixel weighting in case of non-uniform 
sky coverage as well as noise subtraction, provided the noise 
statistical properties (correlation functions) are known. It has 
the additional advantages of being simple and rather fast: 
the required computing time practically corresponds to that 
needed to calculate the two-point correlation functions, i.e. 
about 35 minutes on our platform. The number of operations 
is (9(A^pjj^), while optimal methods to measure power spectra 
require matrix inversions tipically implying a number of 
operations 0(N^- ) (Wright 1996). The implementation of a 
new algorithm £r computing N-point correlation functions is 
expected to further reduce the computing time (Moore et al. 
2000). 



3.2. Destriping signal-only maps 

To check for the presence of spurious correlations introduced 
by our algorithm and depending on the underlying signal, as a 
first test we run the destriping code over CMB and Synchrotron 
maps containing no noise. The test confirms our algorithm sub- 
tracts the average signal from the input map when dealing with 
scalar quantities (CMB temperature in our case), whereas noth- 
ing is lost in our Q and U maps. We then analyse the difference 
between the destriped and the original maps as a function of the 
number of iterations of the destriping cycle. Should no resid- 
ual effects be present, these differences would asimptotically 
become null for Q and U maps, and a constant number cor- 
responding to the average of the input signal for scalar maps. 
Indeed, the maximum pixel-to-pixel diff'erence, which depends 
on the level of the underlying signal, decreases with increas- 
ing number of iterations. In particular, the spurious noise in- 
troduced when destriping maps of synchrotron emission at 
22 GHz, i.e. the SPOrt frequency where the Galactic contri- 
bution is expected to be higher, becomes negligible also for 
CMBP after about 20 loops, reaching the tiny level of few nK 
after 30 loops. More details are shown in Table 1, both in terms 
of peak-to-peak amplitude and rms of spurious noise. For com- 
parison, the peak-to-peak amplitude and rms of the underlying 
signal maps are shown as well. 

The same conclusions can be derived by inspecting the an- 
gular correlation functions C^(ff) (X = Q, U) of the spurious 
noise: they are not flat, but deviations from a flat null function 
decrease with increasing number of iteration cycles and can 
be made negligible after about 10 (20) loops when destriping 
CMB (synchrotron) Q, U maps. 

3.3. Destriping noise maps 

The effectiveness of our destriping technique is then tested on 
simulated maps containing only noise, in the following cases: 

- white noise 

- 1 // noise 

- both white and 1 // noise 
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Table 1. Peak-to-peak amplitude and rms of the spurious noise 
introduced when destriping signal-only maps (CMB and syn- 
chrotron), as functions of the iteration loop number, for Q and 
U maps. The last line shows the peak-to-peak amplitude and 
rms of the underlying Q & U signals. 
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Fig. 3. Noise simulated maps before (top) and after (bottom) 
destriping, for the case = 1 .8 ■ lO-'*, in Kelvin. The HEALPix 
parameter A^i/iie is 32. 



The good performances of the technique are made evident in 
Fig. |3] where we show a simulated noise map (both l/f and 
white) before and after destriping. 

One possible way to quantify the residual correlated noise 
after destriping is measuring the fractional excess pixel noise 
with respect to the case of purely white noise. Results are 
shown in Table 2 for two different values of the knee fre- 
quency, corresponding to the SPOrt goal knee frequency, /k = 
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Fig. 6. Increment, due to the presence of residual correlated 
noise, in the rms of measured noise power spectra, after de- 
striping, as a function of the multipole, in percentage of the 
rms of purely white-noise power spectra. 



1.8x10-^ Hz, and the SPOrt orbit frequency, /o = 1.8x10 '* Hz, 
the latter representing a very conservative case. 

Table 2. Excess rms noise due to low-frequency contributions, 
with respect to the white noise level, for pixels of ^ 7°. 
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After Destriping 


1. 


8 X 10-^ 


310% 


6% 
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8 X 10-5 


35% 


< 1% 



Another way to quantify the residual correlated noise is 
measuring and inspecting the two-point correlation functions 
C^{ff} of simulated Q and U noise maps. Averages and Icr 
bands of 500 correlation functions C^{9) measured from maps 
containing both white and l/f noise, before and after destrip- 
ing, are compared to the purely white noise case in Fig.l^for the 
same knee frequencies as in the previous test. The correlation 
functions C^(9) are similar and are not shown. As expected, 
the correlated noise is strongly reduced by the destriping pro- 
cedure, the residuals falling within the statistical error of the 
white noise case for = /if""'. The correlation functions for 
the case of purely l/f noise are similar to the white plus l/f 
case everywhere but at distance, and are not shown. 

If the noise statistical properties are known, its two-point 
correlation function after destriping can be calculated with 
Monte Carlo techniques and subtracted from the measured val- 
ues before integrating, the recovering of signal power spec- 
tra presenting no huge problems. However, the residual corre- 
lated noise also increases the expected error on measured quan- 
tities. In our case the rms of correlation functions measured 
from many noise maps is increased by roughly 4% (11%), for 
/k = 1.8 X 10-^ Hz (/k = 1.8 X 10-"* Hz), with respect to the 
case of purely white noise. 



For completeness, and to ease the comparison with other 
methods, average noise power spectra, obtained from the cor- 
relation functions used to make Fig. |3 are shown in Fig. |5j 
though they carry the same information as the correlation func- 
tions. As expected, even for the largest knee frequency the Cf 
power spectrum of the residual noise is close to that of purely 
white noise. The region of low multipoles is the most sensitive 
to low frequency residuals, some contributions being always 
found here also after the application of other destriping tech- 
niques (Maino et al. 1999). The excess noise, in percentage of 
the measured white noise level, is shown in the bottom panel 
of Fig.|5las a function of the multipole. The theoretical white 
noise level can be calculated as follows: 

Cr = (24) 

where {cr^} represents the average noise variance per pixel. Our 
measured value is equivalent to its expectation of 0.37 yuK^, 
again confirming that our map-making algorithm is correct. 

Finally, Fig.|6lshows the increment, due to the presence of 
residual correlated noise, in the rms of measured noise power 
spectra, after destriping, in percentage of the rms of purely 
white noise power spectra, again as a function of the multipole. 

4. Conclusions 

We presented a new iterative destriping technique working di- 
rectely on Q and U Stokes parameter data as well as on scalar 
quantities. As for other methods, in order for the technique to 
work properly the receiver knee frequency must be lower than 
the experiment modulation frequency, and there must be suf- 
ficient overlap between different orbits. On the other hand, no 
requirements on the statistical properties of the noise are neces- 
sary here. Moreover, our observing strategy preserves average 
Q and U values in the destriped maps, though the map average 
is still lost when dealing with scalar quantities. 

The performances of the technique, studied on simulated 
SPOrt data by analysing both the measured C^'^(6) two-point 
correlation functions and the measured Cf^ power spectra, are 
comparable to those of other methods based on minimisa- 
tion. Computer time is not an issue, no matrix inversions being 
involved. 

Power spectra are recovered from the measured two-point 
correlation functions by integration, this method being imple- 
mented here for the first time for polarization data. It has the 
advantages of avoiding edge problems arising when using the 
ANAFAST code and to allow for pixel weighting in case of non 
uniform sky coverage. 
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